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The neutrino mixing phenomenon, was firstly studied in the context of quantum mechanics [U 
and more rece ntly analyzed in the framework of the quantum field theory (QFT) formalism [l( 

(Ni : [il,[i3,[ii,[iS[2al2i|. 

The recent experimental achievements proving neutrino oscillations [2 3 . |23[ and the progresses in the QFT under- 
standing [2^ [2l| of the neutrino mixing seems to indicate a promising path beyond the Standard Model of electro- weak 
^ interaction for elementary particles and a possible link between high energy physics and cosmology (2^ . [25| . In this 
^ , paper, indeed, we show that the energy content of the neutrino mixing vacuum condensate [2^ can be interpreted 



^ ' as dynamically evolving dark energy [25|| that, at present epoch, assumes the behavior and the value of the observed 

f~| , cosmological constant. We compute such a value and show that, above a threshold, it is slowly diverging and its 
derivative with respect to the cut-off value goes actually to zero (cfr. Fig. (2)), which allows to use the cut-off at its 

I ' Planck scale value. 

^ I Our result links together dark energy with the sub-eV neutrino mass scale. The link comes from the neutrino- 

^\ ' antineutrino pair vacuum condensate due to the mixing phenomenon. 

T-H , This fact is crucial from a genuine experimental point of view since, up to now, none of the exotic candidates 

' for dark matter and dark energy, has been detected at a fundamental level. Considering neutrino mixing vacuum 

\ condensate as the source of dark energy fits with the conservative view by which only actually observed ingredients 

' as gravity, radiation, neutrinos and baryons are taken into account. 

\ The layout of the paper is the following. In Section II, we outline the neutrino mixing formalism in Quantum Field 

r Theory. In Section III we compute the neutrino mixing contribution to the dark energy in the case of two generations. 

' The case of the three fiavor fermion mixing is analyzed in Section IV and conclusions are drawn in Section V. 



II. NEUTRINO MIXING IN QUANTUM FIELD THEORY 

For the reader convenience the main features of the QFT formalism for the neutrino mixing are here summarized. 
For the sake of simplicity, we restrict to the case of two fiavors. Extension to three flavors ^] can be also considered 
(for a detailed review see [T9|). 

The Pontecorvo mixing transformations for two Dirac neutrino fields are 

Ve{x) — i'i{x) cob9 + V2{x) smO 

'^fiix) ~ — (^1(0;) sinfl + i^2(a;) cos6' , (1) 

where Veix) and ^^{x) are the fields with definite flavors, is the mixing angle and vi and V2 are the fields with 
definite masses mi ^ •m2'. 

v,{x) = -^Y. <M + "-k,,: * = 1, 2, (2) 

k,r 

with al^^{t) = i e-''^'^.'*, /J^^^t) = /^^^^ e*"*" '*, and ujk,, = ^/k^ + ml The operators a^. ^ and i = 1, 2 , r = 1, 2 
annihilate the vacuum state |0)i,2 = |0)i |0)2: Q^k J0)i2 = /3ki;|0)i2 — 0. The anticommutation relations are: 
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[iyrix),iyfiy)}^^^^ = 53(x-y)5„0%, with a,/3 



.4, and 



with i,j — 1,2. All other anticommutators are zero. The orthonormality and completeness relations are: u^^jWk i ~ 



The mixing transformation Eqs.|T]) can be written as [lO|: 

,,:ix) = Gg\t) u'^ix) Ge{t) 
K{^)=G,\t)v^{x) Gs{t) 



(3) 



where the mixing generator Ge{t) is given by 



Ge{t) — exp 



/ d^x (^iyl{x)i^2{x) - iyl{x)iyi{x) 



(4) 



At finite volume, Gg{t) is an unitary operator, Gg ^{t) = G^g{t) — Gg{t), preserving the canonical anticommutation 
relations; Gg^{t) maps the Hilbert spaces for free fields Tii.2 to the Hilbert spaces for interacting fields TCe,^i'- Gg^{t) : 
Til, 2 ^ 'He.fj.- In particular, for the vacuum |0)i,2 we have, at finite volume V: 

|0(i)>e,p=G,-'(t) |0)l,2. (5) 

|0)e.p is the vacuum for Tig./j, which we will refer to as the flavor vacuum. In the infinite volume limit the flavor 
vacuum \0{t))e,ii turns out to be unitary inequivalent to the vacuum for the massive neutrinos |0)i,2 |10|. This can 
be proved for any number of generations [l5|. The non-perturbative nature of the flavored vacuum for the mixed 
neutrinos is thus revealed. 

Due to the linearity of Gg{t), we can define the flavor annihilators, relative to the fields Veix) and Vfj,{x) at each 
time expressed as (we use (cr, i) = (e, 1), (/i, 2)): 



^ G~,\t) alAt) Goit), 
PLAt) ^ G-,\t) PIS) Gg{t). 
The flavor fields can be expanded in the same bases as v^: 



Va{yi,t) = 



1 



k,r 



(6) 



(7) 



The flavor annihilation operators in the reference frame such that k = (0, 0, |k|) are: 

a^M = cosflaLiW + {\U^\ al^t) + e^Vil /3:L2W) 

«k,^W = cos^ai;^^^ - sine {\U^\ a^^it) - ^ ^X^)) 

P-uAt) = cos^CkaW + sin^? (|t/k| /3!:k,2W - \Vu\ all,{t)) 

Ck,pW - cos0/3:k,2W - sin^? {\U^\P-i.At) + e'' l^k| 



(8) 



with e*" = (—1)'^ and 

|[/k| = Uk^.wLj 
with i,j ~ 1,2 and i ^ j. We have: 



-k,i'^-kj ) 



|Vk| 



rt r 
^k,l'^-k,2 



k,2'"-k,l 



(9) 



|C/kh 
l^kl 



mi 



\ 2uJkS 
' LOkA + mi 



2uj 



i^k,2 + m2 

20Jk,2 

^fc,2 + m2 

2uJk.2 



{uJk.i + mi){ujk,2 + "12) 

k k 



(Wfc,2+"l2) {uJk,i+mi) 



(10) 



|(7kP + ll^kP = 1. 



(11) 



Figure 1; The fermion condensation density |V(p, a)|^ as a function of p for a = 0.98 (solid line) and a = 0.5 (dashed hne). 



The condensation density is given by 

^L>lc,.|0)e.p = eAf^lPV^PlMe,^. = sin^ |VkP , Z = 1, 2 . (12) 



The Bogohubov coefficient |Vk| appearing in the condensation density can be written as a function of the dimen- 
sionless momentum p = , and dimensionless parameter a = LuimHh)— as foUows, 



\V{p,a)\^ = - \l^—^^===\ . (13) 




From Fig.(l) we see that the effect is maximal when p = 1, and \V\'^ goes to zero for large momenta (i.e. for 
|k|2 » mim2 ) as \V\^ « 

Since the experimentally observed neutrinos are always extremely relativistic, the value of \V\'^ is very small. Only 
for extremely low energies (like those in neutrino cosmological background) might be large and account for few 
percent. 

In the next Section we will show that the mixing of neutrinos may contribute to the value of the dark energy exactly 
because of the non-zero value of |Vkp: its behavior at very high momenta, together with the Lorentz invariance of 
the vacuum condensate at the present time, can be responsible of the very tiny value of the cosmological constant. 



III. NEUTRINO MIXING AND DARK ENERGY 



Experimental data indicate that the today observed universe can be described as an accelerating Hubble fluid where 
the contribution of dark energy component to the total matter-energy density is Q.\ ~ 0.7 (see the discussion in Section 
V). Moreover, the cosmic flow is "today" accelerating while it was not so at intermediate redshift z (e.g. 1 < z < 10) 
where large scale structures have supposed to be clustered. Thus, physically motivated cosmological models should 
undergo, at least, three phases: an early accelerated inflationary phase, an intermediate standard matter dominated 
(decelerated) phase and a final, today observed, dark energy dominated (accelerated) phase. This means that we 
have to take into account some form of dark energy which evolves from early epochs inducing the today observed 
acceleration. 

In this Section we show that the energy density due to neutrino mixing vacuum condensate can be interpreted as 
an evolving dark energy that at present epoch has a behavior and a value compatible with the observed cosmological 
constant. 

The calculation here presented is performed in a Minkowski space-time but it can be easily extended to curved 
space-times. When particle mixing and oscillations in curved background is analyzed, neutrino mixing, and in general 
particle mixing, gives a time dependent dark energy leading, however, to the same final result: the today observed 
cosmological dark energy value can be recovered. 

Let us calculate the contribution of the neutrino mixing to the vacuum energy density. 

As well known (26j , the Lorentz invariance of the vacuum implies that the vacuum energy-momentum tensor is equal 
to zero: T^^^'^ = (0| : : |0) = 0, (as usual normal ordering is denoted by the colon : ... :). The (0,0) component of 
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the energy-momentum tensor density Tqo{x) for the fields vi and i/2 is then 

: %oix) := ^ : (*m(x)7o do *m(a;)) : (14) 

where ^'m = (2^1,1^2)- In terms of the annihilation and creation operators of fields i^i and 1^2, the (0,0) component of 
the energy-momentum tensor Tqq = / (PxToo{x) is given by 

■ ^(°) ■-T.f '^'kt... (<,ak, + PX^PL^,) , (15) 

with i = 1,2. Note that is time independent. 

In the early universe epochs, when the Lorentz invariance of the vacuum condensate is broken, p™" presents also 
space-time dependent condensate contributions. This implies that the contribution of the neutrino mixing to 
the vacuum energy density is given by computing the expectation value of T^^ in the flavor vacuum \0{t)) ^ 

pzi = ^ %o eAm\ E ■■ ■■ low)^./^ ■ (16) 

i 

Within the QFT formahsm for neutrino mixing, we have 

e,^(0| : T°° : |0)e.^ = eAm\ ■■ T^°) ■■ |0(i))e,M (17) 

for any t. We then obtain 



mix 
rvac 



E / 7^^M(e,,(0|<X.|0)e,,+e.,(0|/3;;XjO)e,,), 



and then 



^ 4sin2 e / ^ (cfe,! + Loua) m\ (18) 



which, introducing the cut-off K, becomes 

2 

pZc ^-sin^e dkk^{ujk,i+LUk.2)\Vk\^. (19) 
Jo 

In a similar way, the contribution p""^ of the neutrino mixing to the vacuum pressure is given by the expectation 
value of T^^^ (where no summation on the index j is intended) on the flavor vacuum |0)e,;i: 



mix 

~ ~v 



Pvac 



:^'toe,^(0|E:?^S:|0)e,^ . (20) 
Being 

: E / ^ + P%P-K.) > (21) 



in the case of the isotropy of the momenta we have T^^ = T^^ = T'^^, then 

2 . 



O TT 





1 


1 
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(22) 



From Eqs. pQl) and (|22p we have that the adiabatic index is w = P^ac I Pvac — 1/3 when the cut-off is chosen to be 

K 3> mi, 7712. 

The values of p""^ and p™^ which we obtain are time-independent since we are taking into account the Minkowski 
metric. Considering a curved space-time, time-dependence has to be taken into account but the essence of the result 
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™" comes from 



is the same. At the present epoch, the breaking of the Lorentz invariance is neghgible and then p; 
space-time independent condensate contributions (i.e. the contributions carrying a non- vanishing ^ — (LUk,kj) 
are missing). That is, the energy- momentum density tensor of the vacuum condensate is given by 



(2^ 



,3 e,P\ 



(23) 



Since rj^^ — (img(l, — 1, — 1, — 1) and, in a homogeneous and isotropic universe, the energy-momentum tensor is 
T^i, = diag{p ,p ,p ,p), then, consistently with Lorentz invariance, the state equation is p™" 



Pa 



This means 



that the vacuum condensate, coming from neutrino mixing, contributes today to the dynamics of the universe by a 
cosmological constant behavior [2^ . ExpUcitly, we have 



Pa 



2 . 2 
— sm 

77 



Solving the integral, we obtain 

?| ('77-2 ^ m\)k 

( 



mix 

Pa 



2 . 2 
— sm 

TT 



fc2 



^m,2\/k? + m\ 



K 



dkk^ 



2(r7i2 - mi) 



(2m^ — 2mfm2 + m\m 



m\) log 



log k 



miy'k^ 



rrin 



(24) 



(25) 



that is 



Pa 



— sm 

TT 



2(m2 - TOi) 
■\/ — 



TO]^ arctan | 



-.K 



1712 arctan 



-.K 



{2m\ 



2 2 
ra^ m 



(2m 



2m2mi 



m\m\ — m\) log ( X 



l"-2 - "12) log (^^■'^ 

— (2m^ — 2ni\m2 



nilnii, — TO2) log (7712) 



(2to4 



2m\mi + m\m\ — ) log (toi) |. 



(26) 



The plot of p™*^ as function of the momentum cut-off K (Fig. 2) shows that for K at Planck scale a value of p™' 
is obtained, which is in agreement with the observed value of cosmological constant. 

To better understand the meaning of Eg. ([26|) . we report the behavior of p™" for K ^ rni, 7712: 



Pa 



arctan 



V mi / 



mi 



m,2 arctan 



\ 1712 ^ - 



{2mf — 2m\m2 + m\in\ — m\) log (m2) — (2m,2 



2r7i\mi + rn^m"^ — m\) log (mi) 



+ {7n{ + 7n\ + 2m\ml - 2mlm2 - 2m^mi) log {2K) |. 



(27) 



This shows that the integral diverges in K as rTij log(iir). As shown in Fig. 2 the divergence in K is smoothed by 
the factor mj. For neutrino masses of order of 10"^ eV we have p™'^ = 5.4 x 10~^^GeT^'* for a value of the cut-off 

of order of the Planck scale K — lO^^GeV^. From Eq. ipT)) one also sees that ''"''^j^j^^^ oc — > for large K. An 



interesting question to ask is how the result p™*^ oc m* log(ii'), directly obtained in our approach, is related to the 
conjecture y!f\ that the small value of the cosmological constant pA oc (lO^^ey)^ is associated with the vacuum in a 
theory which has a fundamental mass scale m ~ lQ~^eV . 

We observe that, since, at present epoch, the characteristic oscillation length of the neutrino is much smaller than 
the universe curvature radius, the mixing treatment in the flat space-time, in such an epoch, is a good approximation 
of that in FRW space-time. The central result of this paper is: the vacuum condensate from neutrino mixing can give 
rise to the observed value of the cosmological constant. Exotic components to dark energy are not necessary in this 
approach. 
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Figure 2: The neutrino mixing dark energy as a function of cut-off K. 



IV. THREE FLAVOR FERMION MIXING 



The above result can be easily extended to the three flavor fermion mixing case. The Lagrangian density describing 
three Dirac fields with a mixed mass term is: 



(28) 



where ^'J — [v^, t'p, J^r) and M = is the mixed mass matrix. Among the various possible parameterizations of the 
mixing matrix for three fields, we work with CKM matrix of the form: 



C12C13 S12C13 Sise 



= W *m(a;) -S12C23 - ci2S23Si3e* C12C23 - si2S23Si3e' S23C13 *m(a;) 



S12S23 - Cl2C23Sl3e** -C12S23 - Sl2C23Sl3e"' C23C13 



(29) 



with Cij — cos 9ij and Sy — sin 9ij, being 9ij the mixing angle between i/i, Vj and '^J-,-^ — {vi,V2tVz)- Using Eq. (P5)) . we 
diagonahze the quadratic form of Eg. ([28)1 . which then reduces to the Lagrangian for three Dirac fields, with masses 
TOi, m2 and TO3: 



(30) 



where M^j = diag (mi, 1112,7113). 

The mixing transformation can be written as v'^{x) = Gg^ {t) vf-[x) Geit)^ where (cr, j) = (e, 1), (/x, 2), (r, 3), and 
the generator is now 



where 



G0(t) =G23(i)Gl3(i)Gl2(i): 



12 / (f X {v\(x)v2{x) - vl{x)vi{x) 



Gi2{t) = exp 

G23(i) = exp 6*23 j d^x {h'l{x)v3{x) ~ vl{x)v2{x)^ 
Gi3{t) = exp [(?i3 f d^x (^i^l{x)i^3{x)e-'' ~ 4{x),yi{x)e'' 



(31) 

(32) 
(33) 
(34) 



The vacuum for the mass eigenstates is denoted by |0)„i and the flavor vacuum is given by \0{t)) f = Gg ^{t) |0), 
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The flavor annihUation operators in the reference frame k = (0, 0, |k|) are: 

aL,ei*) = C12C13 + S12C13 {p^l a^^it) + e^\V^^\ /3l,2W) + e-'' (iC/fal ^.sW + el^^ial /^l.sW) (35) 

= (C12C23 - e''' S12S23S13) aL,2W - (si2C23 + e'^ C12S23S13) (it^J'al aL,i(i) " £^^121 /3-k,i(i)) 

+ a23Cl3 (lC/2'^31 "k,3(*) + en^23l /3l,3 W) ' (36) 
«lc,r(^) = C23C13 QL3(i) " (ci2S23 + e^' .S12C23.S13) {\U^s\ ^1,(1) - e'-|1^23l /?l,2(i)) 

+ (S12S23 - e^' C12C23S13) (\U^3\ <,i{t) - e^\V,^,\ PX,{t)) , (37) 

/3^k,e(i) = C12C13 /3!:k,i(i) + ^i2Ci3 (|C/f2l /3:k,2(i) - e'\Vh\ all,{t)) + e^' S13 {\U^s\ dL^^^sit) - e^\V,^,\ all,{t))m 

P-Ki^it) = (ci2C23 - e-'' S12S23S13) /3:k,2(i) - (si2C23 + e"'* C12S23S13) (|C/i''2l /?-k,i(*) + \Vi2\ <i(i)) + 

+ S23C13 {\U^s\ PLu,3{t) - e'- 1^2^31 <3(i)) - (39) 

P-Krit) = C23C13 /?!:k,3 " (^12^23 + e"'^ S12C23S13) (|C/2^3l /?-k,2(t) + £11^231 <2(*)) 

+ (S12S23 - e-'* C12C23S13) (|C/^3l /?!:k,i(i) + ^11^3 1 <i(i)) • (40) 

These operators satisfy canonical (anti)commutation relations at equal times. and are Bogoliubov coeSicients 
defined as: 



l^k, 1 + [::] (41) 

\U^/ + \V^\' = 1 (43) 
where i,j = 1, 2, 3 and j > i. The condensation densities are different for particles of different masses: 

Af^ = f{0{t)\N^mt))f = f{m\N^n^(*))f = |-«i2C23 + e'*ci2S23ai3r IV^^l' + slscl, \V^,\\ (45) 

A/-^ = /(0(i)|JV^f|0(i))/ = /(0(t)|Ar^'n0(t))/ = |-Ci2S23 + e'*Si2C23Sl3|^V^23l'+|si2S23 + e'^Ci2C23Sl3n 

In this case, at present epoch, the contribution given to the dark energy by the neutrino mixing is 

p-^- = 2 r^^^2ri!i(^2^^2^|^k|2^^2^|^k|2)^i^/|_^^^^^^ 

(|-Cl2S23 + e''' 5120238131^^231^+ |si2S23 + e'''ci2C23Sl3nV"i^|^) , (47) 



_|_ "^3 

which can be written as 



PT^ = ^ r*fc^{^(s?24|V^l'2l' + «?3|l^l'3n + ^[(«f2ci3 + C?2«i3«?3) \y^2f ^ s\A3\y^3?\ 

2 

l^~t~h^ r/2 2 2 2 2\i k |2 ^2 2 22 2\i k 1 2 1 l 

+ " [(^12^23 +•^12^23^13) r23l + (•^12*23 + Cl2C23'5l3) l^^sl J > 

- -S12C23C12S23S13C^ ^ dkkH^ \V^^\^ + ^ [\V^,\^ - \V^,\^] | , (48) 
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where = cos 6. We note that py^*" is also depending on the CP violating phase S. Like in the case of two flavor 
neutrino mixing, the integral diverges in K as mf log {K). A value of compatible with the upper bound on the 

cosmological dark energy, is obtained for neutrino masses of order of 10~^eV so the result is essentially the same of 
the two flavor case. 



V. CONCLUSIONS AND DISCUSSION 



The vacuum condensate generated by neutrino mixing can be interpreted as an evolving dark energy that, at present 
epoch, behaves as the cosmological constant, giving rise to its observed value. The result is naturally achieved even 
when a cut-off K of the order of Planck scale is considered. It is easily recovered also for three flavor fermion mixing. 
Such a result links together dark energy with the neutrino masses. Introducing auxiliary flelds or mechanisms is not 
required in our approach. 

A short summary of the observational status of art can aid to clarify the frame for our considerations and results. 
An increasing bulk of data have been accumulated in the last few years. They have paved the way to the emergence 
of a new standard cosmological model usually referred to as the concordance model. The Hubble diagram of Type 
la Supernovae (SNela), measured by both the Supernova Cosmology Project [1^ and the High-z Team [2^ up to 
redshift z ~ 1, was the flrst evidence that the universe is undergoing a phase of accelerated expansion. On the other 
hand, balloon born experiments, such as BOOMERanG [13] and MAXIMA [3l|, determined the location of the first 
and second peak in the anisotropy spectrum of cosmic microwave background radiation (CMBR) pointing out that 
the geometry of the universe is spatially flat. If combined with constraints coming from galaxy clusters on the matter 
density parameter ^Im, these data indicate that the universe is dominated by a non-clustered fluid with negative 
pressure, generically dubbed dark energy, which is able to drive the accelerated expansion. This picture has been 
further strengthened by the more precise measurements of the CMBR spectrurn, due to the WMAP experiment , 
and by the extension of the SNela Hubble diagram to redshifts higher than 1 [3^ . Several models trying to explain 
this phenomenon have been presented; the simplest explanation is claiming for the well known cosmological constant 
A Although the best fit to most of the available astrophysical data (s^l, the ACDM model fails in explaining 

why the inferred value of A is so tiny (120 orders of magnitude lower) compared to the typical vacuum energy values 
predicted by particle physics and why its energy density is today comparable to the matter density (the so called 
coincidence problem). As a tentative solution, many authors have replaced the cosrnological constant with a scalar 
field rolling down its potential and giving rise to models referred to as quintessence [35|. Even if successful in fitting 
the data, the quintessence approach to dark energy is still plagued by the coincidence problem since the dark energy 
and matter densities evolve differently and reach comparable values for a very limited portion of the universe evolution 
coinciding at present era. In this case, the coincidence problem is replaced with a fine-tuning problem. Moreover, it 
is not clear where this scalar field originates from, thus leaving a great uncertainty on the choice of the scalar field 
potential. 

The subtle and elusive nature of dark energy has led to look for completely different scenarios able to give a 
quintessential behavior without the need of exotic components. To this aim, we observe that the acceleration of 
the universe only claims for a negative pressure dominant component, but does not tell anything about the nature 
and the number of cosmic fluids filling the universe [ssj. This consideration suggests that it could be possible to 
explain the accelerated expansion by introducing a single cosmic fluid with an equation of state causing it to act like 
dark matter at high densities (giving rise to clustered structures) and dark energy at low densities (then giving rise to 
accelerated behavior of cosmic fluid). An attractive feature of these models, usually referred to as Unified Dark Energy 
(UDE) or Unified Dark Matter (UDM) models, is that such an approach naturally solves, at least phenomenologically, 
the coincidence problem. Some interesting examples are the generalized Chaplygin gas l37i] . the tachyon field [38| 
and the condensate cosmology [s^. A different class of UDE models has been proposed where a single fluid is 
considered whose energy density scales with the redshift in such a way that the radiation dominated era, the matter 
dominated era and the accelerating phase can be naturally achieved. Actually, there is still a different way to face 
the problem of cosmic acceleration. It is possible that the observed acceleration is not the manifestation of another 
ingredient in the cosmic pie, but rather the first signal of a breakdown of our understanding of the laws of gravitation 
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42^1 . Examples of models comprising only the standard matter are provided by the Cardassian expansion [43] , the 
DGP gravity ,44] . higher order gravity actions |45j, non - vanishing torsion field [46.] , higher-order curvature invariants 
included in the gravity Lagrangian [43], etc.. 

This abundance of models is from one hand the signal of the fact that we have a limited number of cosmological 
tests to discriminate among rival theories, and from the other hand, that a urgent degeneracy problem has to be faced. 
The evidences of neutrino oscillations p^ . [23| and the fact that the vacuum condensate originated by neutrino mixing 
provides contributions to the dark energy compatible with today expected value, as shown in the present paper, could 
contribute towards a solution of such a problem from both experimental and theoretical viewpoints. 
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